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I. INTRODUCTION 



Rashba precession 1 of spins and its manipulation has 
been studied intensively 2 in recent years, aiming to con- 
trol the coherent propagation of electron spins. One goal 
is to realize a spin transistor 3 ; another interesting op- 
tion could be to switch between singlet and triplet en- 
tangled states for directing, for example, noise statistics 4 
or for quantum computing. Despite of major efforts, 
manipulation of coherent Rashba precession could not 
be demonstrated experimentally yet in the polarizer- 
analyzer type transport arrangement. Rashba spin split- 
ting occurs near surfaces in the presence of internal or 
externally applied symmetry breaking electric fields per- 
pendicular to the transport direction 1 . Spins precess 
when injected out of the spin-orbit eigendirections. Any 
dephasing along the structure limits successful transis- 
tor operation. One important dephasing mechanism in 
spin-orbit active structures of more than one dimension 
arises due to momentum randomizing scattering events 
by impurities (elastically) , by phonons (inelastically) 5 , or 
by electron-electron scattering events 6 . One-dimensional 
structures confine the direction of propagation, thus re- 
ducing this source of dephasing mechanism 7 . Already 
in their original proposal 3 Datta and Das therefore sug- 
gested to use clean one-dimensional structures for the 
spin-orbit active medium. The current work focuses on 
quantum wires and assumes absence of momentum ran- 
domizing scattering events. No spin relaxing mechanisms 
(in the sense that off-diagonal entries of the spin density 
matrix decay) will be considered here. 

Most quantum wires of current experiments accom- 
modate more than one transport channel at the Fermi 
energy. This holds true also for single walled carbon 
nanotubes (NT) where at least two flavor channels carry 
(spin) current. In NT spin transport has been estab- 
lished experimentally 8 and studied theoretically 9 . Other 
experimental multichannel systems consist of arrays of 
quantum wires fabricated artificially in parallel at close 
proximity 10 . In any multichannel wire the interesting 
question arises whether Rashba spin phases increase by 
equal amounts along different channels or not. Provided 
the kinetic energy dispersion is strictly parabolic, de- 
scribed by a common effective carrier mass in each chan- 
nel as in most semiconducting quantum wires, one would 



expect equal spin phases in all channels as a result of 
the linear splitting by the Rashba energy 3 , so that their 
probabilities to enter the spin selective drain contact all 
add up. We shall demonstrate that Coulomb interactions 
between charged electrons affect their spin propagation 
properties and ultimately destroy the phase relationship 
of Rashba precessing spins between channels in quantum 
wires (QW) fabricated on the basis of semiconducting 
material. This limits operation of the Datta Das spin 
transistor. In metallic single walled nanotubes, on the 
other hand, with their linear kinetic energy dispersion, 
spins dephase between modes already without account- 
ing for interactions; remarkably, in this case we find that 
the Coulomb interaction suppresses this single particle 
dephasing between the two flavor modes, particularly at 
small doping, thus facilitating coherent spin precession 
along NT. 



II. INTERACTING MULTIMODE QUANTUM 
WIRES 



Before addressing the effect of spin-orbit coupling 
let us first discuss the (possibly screened) Coulomb 
interactions within and between transport channels. 
One-dimensional systems are particularly susceptible to 
electron-electron interactions, even when weak 11 . Nei- 
ther in semiconducting QW 12 nor in NT 13,14 interac- 
tions can be disregarded. Contrary to higher dimensions 
they show up, for example, as nonuniversal power laws 
at low energies. The most convenient theoretical frame- 
work is the Bosonic Tomonaga-Luttinger liquid (TL) 15 
developed since the seminal work by Haldane 16 . How- 
ever, a priori, it is not clear in how far the TL-model can 
be applied to multichannel situations. Strictly equivalent 
modes of equal particle densities and interaction matrix 
elements, for example, tend to stabilize a gapped charge 
or spin density, non-TL low energy phase 1 due to the 
appearance of relevant (in the sense of a perturbative 
renormalization group treatment) momentum conserving 
intermode backscattering processes. On the other hand, 
most real systems lack such a strict equivalence of modes. 
In multichannel quantum wires particle densities differ at 
given Fermi energy, and even systems of Ref. 10 without 
fine tuning of densities, or the two flavor modes of NT 
at not exact zero doping cease to be strictly equivalent. 
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Here, we therefore focus on the generic case of unequal 
modes which are expected to stay in the gap-less low en- 
ergy TL phase. 

In this phase Matveev and Glazman (MG) have com- 
puted a power law for the density of states of spin-less 
electrons at the end of the n-th mode 18 : 



v n (w) 



,/3„ 



(1) 



generalizing the single channel case 11 ' 16 ' 19 . MG obtained 
the exponents (3 n within a pure plasmon model by mod- 
eling the charge density fluctuations as coupled har- 
monic strings. Plasmon velocities S£ and the normalized 
eigenmodes 7„£, as obtained from the dynamical matrix, 
determine 18 



h'n 



1 + E hne\ 2 se/v n 



(2) 



v n is the Fermi velocity in the n-th channel. This ap- 
proach merely accounts for long wave length charge prop- 
erties and can shown to be equivalent to the random 
phase approximation (RPA) to plasmon velocities 20 gen- 
eralized to coupled modes 21 . It has been pointed out 22 
that the RPA result deserves improvement at small parti- 
cle densities. Moreover, it disregards exchange processes 
and spin, and tacitly presumes 'super' Galilei invariance 
of every channel individually, as discussed below. Spin 
properties and exchange are known to depend on short 
wave length properties of the interaction 23 . 

For given microscopic interaction the parameters of the 
single channel TL model, and therewith sound velocities 
and exponents, have been obtained from homogeneous 
and static susceptibilities 22 , exploiting exact thermody- 
namic relations 16 ' 15 ' 24 . These susceptibilities, in turn, 
can be computed to high accuracy by standard many- 
body techniques from the underlying Fcrmion model, be- 
yond the RPA or pcrturbative accuracies. This way the 
asymptotic behavior of correlation functions has been 
determined in the ID Hubbard model from the Bcthc 
Ansatz ground state energy Eq 24 . In the Galilei in- 
variant charge sector of QW only the compressibility 

K = j is required (N is the particle number 

and L the system length) to fix the exponent parameter 
K p = In spin sector SU(2) spin rotation in- 

variance enforces K a = 1 (cf., e.g., Refs. 15 and 25). As a 
further important property of one dimension charge and 
spin density wave excitations are expected to separate 
and move at different velocities 16 . Evidence for this has 
been found in recent experiment 10 . Spin velocities of QW 
have been deduced from quantum Monte Carlo magnetic 
susceptibilities 26 . Here we generalize this a priori exact 
thermodynamic approach to coupled channels. Spins can 
be incorporated as separate s = ± modes which allows 
to account even for non SU(2) invariant situations, as it 
arises for example in the presence of a Zeeman- field 27 . 

From the microscopic point of view we consider the ID 
system 
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Fermi annihilation operators c n ^, s refer to wave vector k 
and spin s of mode n. Assuming electron wave functions 
of the product form <~ ip n {x±_)c lkx / '\[L (for a discussion 
of this product assumption cf. Ref. 22) we obtain 

V nin2n3 „ 4 (qr) = Hi- J dx ± J dx' ± x (4) 

from the 3D-Coulomb interaction between electrons, e is 
the dielectric constant of the material surrounding the 
wire. The V nin2n3n4 can be expressed analytically for 
many cross sections of physical relevance 28 ' 29 . Some of 
the V ni n 2 n 3 n 4 vanish by symmetry: for example, angular 
momentum conservation on a cylinder surface of NT re- 
quires ri\ + ri2 = «3 + TI4, or n, must be even in QW 
of mirror symmetric cross section. 'Direct' terms may be 
approximated as V^ n , n {q) « £e*K («) or V^, n , n {q) « 
^§-Io(r|q , |)Ko(r|q , |) at q^ 1 larger than the diameter (QW) 
d or the radius r (NT). At small q both reveal the same 
logarithmic increase which eventually will be screened by 
remote metallic gates. In the above formulae \ v and 
are Bessel functions and q — d 2 q 2 /8. 'Exchange' terms 
between the lowest two modes of parabolically confined 
QW are Vg? 2 (q) = ^ge«(Ki(g) - K (q)) 28 while in NT 
the two lowest degenerate flavor modes have angular mo- 
mentum zero and interact ~ Vmi(?)- 

The Boson model 16 , describing gapless excitations of 
Fermions in a ID wire of length L (periodic boundary 
conditions) in the vicinity of the Fermi energy, takes the 
form 

h = E K } + ^ 2) ) + ik 0^ + ( 5 ) 

when generalized to M modes. Bold letters indicate 
M x M matrices and the components of vectors refer 
to modes 1 < n < M. In Eq. (5) 
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r=± 



describe excitations of right/left (r = +/—) going 
Bosonic density fluctuations g r n{q) at wave numbers 
q 7^ 0. Topological excitations [last term in Eq. (5)] de- 
scribe changes in the ground state energy when particles 
N n or currents J„ are added to mode n at q = 0. This 
second part of Eq. (5) is important in the present context 
and is governed by generalized homogeneous and static 
compressibilities and (Drude) conductivities: 
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and (vj) r 
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respectively. They generalize corresponding TL parame- 
ters for single channels 15 (E is the ground state energy 
of the interacting electron system) and govern the com- 
plete low energy physics. As in the single channel case 16 
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v N and Vj are related to the interactions between density 
fluctuations in Eq. (6), moving in the same 



or opposite 



V {1) /ir= (v N +v,)/2-v 



V^/n = K - v 3 )/2 



(8) 
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directions, assuming e n (k) = e n (— k) and thus Fermi ve- 
locities (v) nn i — v n 5 nn > of equal magnitudes at cither 
Fermi point. 

Importantly, all entries of v N and v., are observable, at 
least in principle, and must therefore agree in Bosonic (5) 
or Fermionic [Eq. (3) together with (7)] representation. 
This allows deduction V^ 1 ^ and by standard many- 
body techniques from the given microscopic interaction 
(4) through Eqs. (7 — 9). Therefore v N and v., (and not 
exponents as often assumed for single channel quantum 
wires) are the principle parameters governing the low en- 
ergy physics in the multichannel case. They establish the 
quantitative link to the microscopic Fcrmion model (3). 
For real quantum wires self-consistent Hartree-Fock 22 ' 30 , 
diagrammatic 23 , or Monte Carlo techniques 26 were used 
to deduce TL parameters from a microscopic model for 
electron-electron interaction potential. As a first ap- 
proach to multichannels we rely on the perturbative ap- 
proximation below, which, when including the Fock term 
(as crucial in spin sector), proves already as superior 
to the often used random phase approximation 20 to TL 
parameters 22 . 

The matrices v N and Vj reflect symmetries of the sys- 
tem under permutations of modes. One important case 
are equivalent modes when (v) nn > = vS nn i in Eq. (6), 
and v N and Vj both are cyclic matrices. Then Eq. (5) 
can be diagonalized and the Bogoliubov transformation 
solving for plasmon velocities and exponents can be car- 
ried out separately in each of the resulting independent 
'normal mode' TL (two equivalent channels, for example, 
can be separated trivially into independent gi + qi and 
Qi — Q2 modes). Another symmetry is Galilei invariance 
(observed in QW and NT) for which Trt>., = J2 n v n 31 
stays independent of interactions. The higher symmetry, 
where the ground state energy changes only by the trivial 
kinetic part, independently of the interaction strength, 
under boosting any individual channel, J n — > J n + 5 n , 
we call 'super' Galilei invariance. It implies v 3 = v, i.e., 
= which in general is not observed by Eq. (3) 
though tacitly assumed often in theoretical work. 

Of central importance is the single particle density ma- 
trix (tp n (x)ipni {0)} which can be evaluated from Eq. (5) 
using Boson operators a nq = y / 2n/L\q\ J2 r ®( r <l)Qrn{q) 



with \a nq , a\, ql ] 



Snn'Sqq', generalizing the single chan- 



nel case 16 ' 15 . New Boson operators b nq eventually di- 



agonalize h[ 1] + H q *> . They are obtained via a Bogoli- 
ubov transformation reading in the multicomponent case 

(vu) ( !)• Here ' the M x M matrices u and 
v must satisfy (w[vuv _1 (wvuv _1 ) t — v(wv) t ]);;' = Sw 
where w = (vuv x u — v 2 ) -1 and t indicates the trans- 
pose. This condition is fulfilled when 



r(2) 



(?) 



provided R is orthogonal, and cu> — 5u> cosh z9; and 
sw = $w sinhi?; arc diagonal matrices, 1 < 1,1' < M. 
The \M{M — 1) real parameters of R and the M angles 
$1 are found from the ^M(M + 1) equations 
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Here, the tilde connotes rotated symmetric matrices, 
A = R UR = A\ 

Finally, the M eigenvalues of 



(c s)/^ + y (1) /^ 
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(12) 



u = R c R 



and 



v = R s R , 



(10) 



with the outcome of Eq. (11) inserted for c, s, 
and R, are the (generalized) plasmon velocities of 

J2 q =io (-^4^ +-^j 2 ^)- The same sequence of transfor- 
mations does not, in general, diagonalize the topological 
excitations which explains why coupled TLs cannot nec- 
essarily be decomposed into independent 'normal mode 
TLs'. 

Equations (7 — 12) allow to calculate the (asymptotic) 
power law behavior of any Fermion function, cf. Eq. (18) 
below. This completes the solution for low energy prop- 
erties in multi-mode electron liquids in ID. In the (super- 
Galilei invariant) special case, (V^)nn' = (V' 2 ')™' = 
Vo for all n and n' , the ensuing open boundary exponent 
as well as the plasmon velocities Eq. (12) agree with the 
results of Ref. 18. As stated above, the present approach 
includes spins through separate modes s = ±1. 

With SU(2) spin rotation invariance charge-spin sepa- 
ration continues to occur in multichannel systems, so that 
introducing the additional index v = p,a for the charge 
and spin sector, respectively, renders v N and v 3 block di- 
agonal. Furthermore, the chiral SU(2)xSU(2) symmetry 
of the fermionic model at low energies provides separate 
spin rotation invariance of right and left movers, giving 
rise to v Na = v Ja in spin sector. Very strong spin-orbit 
coupling will spoil charge-spin separation in general. In 
QW, however, charge and spin are mixed only to the or- 
der C(a 5 ) 32 which therefore is rarely expected to be of 
importance [cf. Eq. (13) below, a is the spin-orbit cou- 
pling parameter] . In NT, on the other hand, charge-spin 
separation is already broken to the order 0(a). Here, 
the very small value of a, expected from the small car- 
bon mass, leaves this breaking less important and justifies 
calculating Rashba precession lengths to the leading or- 
der in a. With charge-spin separation Rashba precession 
is controlled solely by the spin sector v — a where, from 
now on, the index n runs only over M spatial modes. 



III. RASHBA SPIN PRECESSION LENGTH 

We address the physical situations referred to in the 
Introduction by studying two model systems: QW for 
quantum wires and NT for nanotubes. In QW the sin- 
gle particle kinetic energy dispersion may be regarded as 
e n W (^) = £ n + k 2 /2m e where e„ are subband energies and 
to is the effective mass. The Rashba spin-orbit energy 



H so = a(a x k z - a z k) , (13) 

deemed as independent of the mode index, adds a lin- 
ear contribution ±ak and splits e^ w (fc) into two equal 
parabolas intersecting at the origin (this statement holds 
strictly true only for not unrealistically strong spin or- 
bit coupling strengths, a <C v n , as discussed above 32 ' 33 ). 
Both branches are indexed by the electron spin projec- 
tion perpendicular to the wire axes, taken as the z- 
direction for convenience (the £-axes points along the 
wire). Within each mode, at the Fermi energy, the 
spin splitting yields a finite difference of Fermi momenta 
fcnT — k n i = 2am c which, as important property of the 
QW case, does not depend on n. Spins initially polar- 
ized along the x-axes therefore precess once over a length 
47r/|fc„| — k n \\ = 27r/am , equal in all channels. 

The other system we consider are metallic carbon nan- 
otubes (NT). Here, spin-orbit coupling 34 arise predomi- 
nantly from the curvature of the tube surface, cf. Ref. 35; 
in flat graphite layers it vanishes by mirror symmetry. 
The kinetic energy dispersion may be taken as e^ T (fc) = 
±{e„ + «F[|fcF — k\Q(k) + \k F + fc|0(— fc)]}, disregarding, 
for simplicity, parabolic parts of the dispersion close to 
the subband bottoms at e n . Here, kp denotes the Fermi 
momentum at zero doping and i>f the Fermi velocity. For 
every n there exist two different flavor branches b = ±1, 
depending on the magnitude |fc| ^ |fcp|, of opposite veloc- 
ities at given sign of k. Adding H so from Eq. (13) yields 
differences k n ^ — fc„j = 2a(e F - e n + bvpk-F)/(v F — a 2 ) 
at the Fermi energy, and therefore Rashba precession 
lengths, depending now on n and, additionally, on b. 

So far we considered independent electrons for QW 
and NT. It has been shown recently that electron- 
electron interactions influence and actually enhance 
Rashba precession in 2D structures 36 and in single trans- 
port channels 32 . To calculate Rashba lengths on the ba- 
sis of the multi channel TL model (5) requires expressing 
Eq. (13) in Bose variables. As in single channels 32 H so is 
found to be proportional to the spin currents J n ,<r- The 
expressions 

C = -«"»e^U»,, (14) 

n 

and 

H so = -av F > J n . a , (15) 

„ v ¥ -a z 

however, differ in the two cases, QW and NT, as can 
be checked in the limit of vanishing interaction, when 
(tw)nn' = {v Ja )nn' = v n S nn > in Eq. (5) reproduces the 
Rashba precession lengths 27r/|fc„| — fc„jj in individual 
QW or NT channels, respectively (H so is a single particle 
operator and does not depend on the Coulomb interac- 
tion). 

With Eqs. (14) or (15) we are now in the position to 
calculate 

fm = iE((ttW+^(i))W T (o)+4(°))) . 

nn' 

(16) 



employing the well known Boson representation for Fermi 
operators ip ns (x) 16 . In the spirit of the Datta-Das setup 3 
the quantity f(L) can be interpreted as the probability 
amplitude to measure an electron at x = L with its ini- 
tial spin polarization parallel to the x-axes, provided it 
was injected at x = with equal probability amplitudes 
into all occupied channels. The result can be represented 
as 

f(L) = w(£)cos(7rL/A„) . (17) 

nn' 

As indicated in the previous section the asymptotic power 

law decay of the two point function g nn '(x) ~ 

can be expressed through the solutions of Eq. (11) for R 

and 

firm' = 5 nn > + 2 RnjRn'j Smh 2 "d j . (18) 

3 

Equation (18) generalizes the known result for the elec- 
tron Green function of single channels 11 ' 16 . By virtue 
of conformal invariance the time dependence can be de- 
duced, yielding the open boundary exponent 

/^-l + ^TRL'C- 2 * 5 "' (19) 

n' 

for the density of states [Eq. (1)]. 

Equation (17) allows to read off the inverse precession 
lengths 

1 am c v— -i x /r .^\ 

^qW = — ZJ?3<T )nn'Vn' (20) 
n ' 

and 

1^ = ^ E (2D 

n 

e n / <|ep| 

e F - e„/ + bv^kp 
v F — a 1 

after which spins reverse their polarization in channel n. 
We see that the A„ are governed by the matrices v 3a 
of spin conductivities. Without electron-electron inter- 
action (ujcr)nn' = v n 5 nn i so that A^ w = A stay equal 
in all channels of QW [Eq. (20)]. This would yield op- 
timum transistor operation as discussed before. Inter- 
actions, however, alter the diagonal entries and, addi- 
tionally, generate off-diagonal entries in v Ja , describing 
the coupling between channels. In general A^ w become 
n-dependent so that different channels dephase. Remark- 
ably, as demonstrated now, interactions reduce spin de- 
phasing in the two lowest flavor modes of single walled 
NT [Eq. (21)], compared to its magnitude in the absence 
of interactions. 



IV. PERTURB ATIVE ESTIMATE 

Sufficiently weak interactions can justify the perturba- 
tive estimate to v Ja . Crucial is the magnitude of nonzero 
Fourier components, which even in metallic NT are small 
compared to the Fermi velocity, since the (large) q = 



A 



component of the interaction does not affect the spin 
sector and is properly accounted for on the RPA level. 
Imposing SU(2) symmetry [naive low order perturbation 
theory violates SU(2) invariance 22 ] results in 

(^N(r)nn' = (v.ja)nn' = [ v n {Vnnnn (2fc n ) 

(k n - kj)}/2ir]6 nn' 

- [V nn >nn>{kn - M/M 1 - 5 nn <) ■ (22) 

Here, k n denotes the Fermi momentum of mode n 37 . 
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FIG. 1. Difference \\ 2 -\i\/\ {0) versus Fermi energy ep /u)o 
in units of the subband energy loq for two channels, cf. main 
text. In QW (a) interactions always cause dephasing between 
both modes while in single walled NT (b) dephasing between 
the two flavor modes can be significantly suppressed at small 
doping compared to the noninteracting case (dashed). 

Figure 1 (a) shows the resulting difference | A2 — Ai | / A' ' 
of Rashba lengths relative to this length in the 
absence of interactions in a two-channel QW of width 
d = cib/V2 (cib is the Bohr radius). This phase differ- 
ence can exceed 30% and decreases only when ef and the 
carrier density increase so that the interaction strength 
diminishes. We see that dephasing never vanishes within 
the regime of validity of the perturbational approach at 
not too small carrier densities in the second subband, 
£f > 1.15co>o [note that for harmonic confining potential 
the third subband (not included here) becomes occupied 
above ep/^o > 2]. Experimentally, it might be possible 
to disentangle two Rashba periods and thus, by monitor- 
ing their dependence on the interaction strength (carrier 
density), verify the predicted dephasing mechanism when 
two channels are occupied. For the Datta Das transistor 
based on semiconductor QWs this result clearly suggests 
to use only the lowest subband. 

In Figure 1(b) the difference of Rashba lengths of 
the two lowest (degenerate) flavor modes of metallic 
single wall arm chair (m, m) NT is seen for m = 5 
[subband energies ei = £2 = 0, the Fermi momenta 
fcf 1 } = vaj\f%r ± e-p/v-p, r is the tube radius, and 

e 2 /evF i=a 2.7 14 ], relative to the splitting A^ at zero dop- 
ing. Dephasing can shown to be suppressed logarithmi- 
cally ~ cf/[ci — C2 ln(eF/wo)] + 0(tp) at small doping (ci 
and c 2 depend on the tube radius and c 2 vanishes with 
the interaction), compared to the interaction free case 
(dashed) where c\ = to/2. This suggests the use of sin- 
gle walled NT close to the neutrality point 38 for coherent 
Rashba precession along both flavor channels 39 . 



V. SUMMARY 



We have calculated Rashba spin precession lengths A„ 
in multimode quantum wires, accounting for electron- 
electron interactions, using the Bosonization method. To 
this end we have developed an in principle exact descrip- 
tion for the TL phase of coupled quantum channels that 
allows the inclusion spin. Generalized charge and spin 
compressibilities and conductivities are identified as the 
key parameters to determine the power law exponents 
and to establish quantitative contact with the underly- 
ing interacting electron model. In semiconducting wires, 
characterized by a parabolic kinetic energy dispersion, 
we find that the A„ become n-dependent, giving rise to 
doubts whether multichannel systems can be used as ac- 
tive part of the Datta-Das transistor. This result could 
explain the up to date lack of successful transistor op- 
eration and clearly suggests to use single channel quan- 
tum wires. In metallic single walled carbon nanotubes, 
on the other hand, we find that dephasing between the 
two flavor modes, arising due to the linear kinetic energy 
dispersion, is suppressed by the electron-electron inter- 
action, particularly at small doping, which could make 
these systems interesting for coherent spin transport. 

I thank Hermann Grabert for many useful conversa- 
tions. This work has been partly supported by the DFG 
under SFB 276. 
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